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The on-shell massless planar double box diagram with an irreducible 
numerator 
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Using a Mellin- Barnes representation, we compute the on-shell massless planar double box Feynman diagram 
with an irreducible scalar product of loop momenta in the numerator. This diagram is needed in calculations of two 
loop corrections to scattering processes of massless particles, together with the double box without numerator 
calculated previously by Smirnov. We verify the poles in e of our result by means of a system of differential 
equations relating the two diagrams, which we present in an explicit form. We verify the finite part with an 
independent numerical check. 



1. Introduction 

We shall consider the following class of Feyn- 
man integrals 
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where the propagator denominators P{ are de- 
fined by P t = k 2 , P 2 = (k+ Pl ) 2 , P 3 = (k + Pl + 

P2) 2 , Pi = (l+pi+p 2 ) 2 , P 5 = (I+P1+P2+P3) 2 , 
P 6 = l 2 , P 7 = (l-k) 2 , P 8 = (k+ Pl +p 2 +p 3 ) 2 , 
d = 4— 2e is the space-time dimension, and the ex- 
ternal momenta pi, P2,P3, and p4 = —pi —p 2 ~P3 
are lightlike. The integrals ([l]) then depend on 
two scales: s = (pi +P2) 2 and t = (p 2 +P3) 2 . We 
shall mainly be concerned with cases where vs is 
either zero or a negative integer. These integrals 
correspond to the diagram shown in figure and 
are important for calculations of two- loop scatter- 
ing amplitudes of massless particles, for example 
inQCD[|,@. 

The scalar double box integral 

If(s,t)= 7(1, l,l,l,l,l,l,0;d), (2) 

was calculated analytically by Smirnov [| as a 
Laurent expansion in e to 0(e ). Subsequently, 
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Figure 1. The planar double box diagram 



Smirnov and Veretin constructed a general algo- 
rithm [ U to deal with tensor double box inte- 
grals. In this algorithm, integrals with arbitrary 
polynomials of loop momenta in the numerator 
are first expressed in terms of scalar integrals 
7(^1, V21 V3, V4, 1/5, v§ , u 7 , 0; d + 2n) with increased 
integer powers of the propagators in shifted space- 
time dimensions J_||- Then, using integration by 
parts identities [ |(| , these scalar integrals are sys- 
tematically reduced to linear combinations of a 
set of master integrals, consisting of I±, a second 
master double box integral, 



Ii{s,t) =7(1, 1,1, 1,1, 1,2,0; d), 



(3) 



and a small number of more simple master inte- 
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grals (see figure ||). Furthermore, they gave an 
explicit analytical expression for the second mas- 
ter integral 7 2 , again as a Laurent expansion in e 
to 0(e°). 



Diag 



Bbox 



L Btri 



Sunset 



jd 

1 Spec 






oo 



Figure 2. Master integrals with less than six prop- 
agators 



It appeared that the algorithm described in 
ref. [ ^| completely solved the problem of calcu- 
lating on-shell double box diagrams. However, as 
was reported by Glover at this workshop [ |j, it 
often happens that in the reduction of a given ten- 
sor integral, the coefficients in front of the master 
integrals I\ and 7 2 are of order 0(1/ e). This is a 
consequence of the fact that in the reduction of 
these integrals it is necessary to reduce the dimen- 
sion down to d = 4 — 2e from at least d = 6 — 2e, 
and in the system of equations for the dimensional 
shift, there are factors of l/(d— 6) sitting in front 
of I\ and Ii. Thus, in order to calculate such 
tensor integrals to 0(e), one would need to know 
1\ and I2 to one order higher in e than they are 
given in refs. [ ||, [3] . 

A typical example is the following integral with 
an irreducible numerator: 



whose reduction to master integrals reads 
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7 7 d rr ( S ,t) =7(1, 1,1, 1,1, 1,1,-1; d), 



(4) 



(6) 

The factors of (d — 4) in the denominators of the 
first two terms on the right hand side of eq. (^) 
are the ones we are concerned about. 

The purpose of this paper is to calculate the in- 
tegral with the irreducible numerator I\ rr directly 
from a Mcllin-Barnes representation. It can then 
be used instead of I2 as a new master integral. We 



3 



check our result in two different ways: firstly, by 
verifying that I± and Ij rr satisfy a system of dif- 
ferential equations, and secondly, by using them 
to compute I± and I2 in d — 6 dimensions, both 
of which are finite, and comparing the result with 
a numerical integration. 

2. Calculation by Mellin-Barnes contour 
integrals 

The analytic structure of the on-shell double 
box is rather simple, since it only depends on two 
scales, and its only thresholds are at s = and 
t = 0. The main difficulty in calculating this 
diagram is to find a way to isolate its infrared 
and collinear divergences. A general recursive 
solution to this problem based on Feynman pa- 
rameters was proposed in ref. [|t]]. However, for 
the analytical calculation we are going to present 
here, it is more convenient to use a Mellin-Barnes 
representation [ ^, || , which enables us to isolate 
the poles in e in a very natural way. 

We derive our Mellin-Barnes representation for 
the two-loop integrals (|l|) by doing the loop inte- 
grations one by one (cf. the off-shell calculation [ 
^)). In terms of Feynman parameters, the /-loop 
can be written as 
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where ^4567 = V\ + ^5 + v§ + v 7 (similar abbrevi- 
ations will be used below), and 

C = —P\x§x 7 — P3X4X7 — P$xc,X7 — SX4XQ . (8) 

By introducing three Mellin-Barnes parameters, 
a, (3 and r, we split the polynomial C into factors: 

ioo ioo 
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x(-sa; 4 a; 6 )' i / 2 - I/1567 - a -' 3 - T r(-a)r(-/3) 
xr(-r)r(i/ 45 67 -d/2 + a + /3 + t). (9) 

After inserting eq. (||) into eq. (|?|) , we evaluate the 
Feynman parameter integrals in terms of gamma 
functions, which gives us the following Mellin- 
Barnes representation for the /-loop: 
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xr(- a )r(-/3)r(-r)r(^567 - d/2 + a + p + T ) 

xT(d/2 - v m7 - a - r)T(d/2 - v i57 - p - t) 
xT(u 5 +T)T(u 7 + a + p + T). (10) 

When this result is inserted into ([!]), the remain- 
ing fc-integral has the form of an on-shell mass- 
less one-loop box diagram with indices v\ — a, 
^2, ^3 — P, v$ — t. We repeat the above steps 
for this fc-integral, using a further Mellin-Barnes 
parameter, a, and finally obtain 

{m d / 2 f (-1)^2345678 \ 
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^ r(a)r(^ 123 8 d/2 a P t (t) 

T(d - fi238 + a + P + t) 
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xT(u 5 + r)r(^ 7 + a + P + r) 
xT(d/2 - i/ 128 + a + r + o)T{v& - t - a) 
xT(d/2-v 23 & + P + T + a)T(is 2 -a). (11) 

In deriving this formula, we have assumed that 
the various parameters are such, that all the ma- 
nipulations we performed are justified. This is 
certainly true if we are able to find a set of straight 
lines, parallel to the imaginary axis, for the inte- 
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gration variables a, /?, a, and r, such that the ar- 
guments of all the gamma functions in it have pos- 
itive real parts. We then define the integrals ([!]) 
for values of the parameters where such contours 
do not exist by analytic continuation. 

Let us now consider the case with the irre- 
ducible numerator, li rT . On the one hand, from 
the definition (|]), 1/5 + v% = 1 — 1 = 0. On the 
other hand, if the real parts of the arguments of 
all gamma functions are positive, then in partic- 
ular i?e(^5 + r), Re{a) and Re(vs — r — a) are 
positive, and therefore i?e(^5 + v%) > 0. Since 
this does not depend on d, it means that in order 
to calculate Ij rr using the Mellin-Barnes repre- 
sentation (|ll|), we must perform an analytic con- 
tinuation not only in d, but also in ^5 or vg. We 
choose vg. Setting i/g = — 1 + 77 and all other z/'s 
equal to one, we get 

(i-K d / 2 ) 2 
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17J.0 1 (— le) 
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- 2e + a 

xT(-l - e - a - r)r(l - r) - e 

xT(-l - e - - r)r(l - 17 - e • 

xr(r; + e - a — — t — a). 

We can make the real parts of the arguments of all 
gamma functions in eq. (12) positive by picking, 
for example, rj — 12y and e = — 12y, where y is 
some positive number much smaller than one, and 
choosing contours for the Mellin-Barnes variables 
defined by: Re (a) = Re{/3) = -y, Re{r) = -1 + 
Ay and Re(cr) — Ay. 

Starting from these values, we first perform an 
analytic continuation in 77 to rj — 0, keeping e 
fixed, and then another one in e to the vicinity of 
e = 0. The procedure for both continuations is 
straightforward: keeping the integration contours 
fixed, we simply have to keep track of the poles of 
the gamma functions, and whenever one of them 



crosses an integration contour, add its residue to 
a list of terms that will contribute to the final 
answer. For example, with the above choice of 
contours, the first crossing happens at 77 = 8y, 
when the pole at t = — 1 + 77 — a crosses the r- 
contour. After taking a residue in one integration 
variable, we continue to follow the poles in the re- 
maining variables, building up a tree of single and 
multiple residue terms. By this procedure, poles 
in e are automatically expressed through singular 
gamma functions multiplying integrals that can 
safely be expanded under the integral sign. 

To C(e°), it turns out that, along with terms 
where there is no integral left, only single and 
two-fold integrals contribute, because terms with 
more integrals are killed by the factor 1 /T(— 2e) in 
eq. ([l2]). In the two- fold integrals, one integration 
can be done by Barnes's first lemma [ fl0| . The 
single integrals that are left can all be calculated 
by closing the contour and summing harmonic se- 
ries. We find the following result: 
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- V - + + 26^ bg(1 + t/g) 

-5(^+7T 2 )l0g 2 (l+i/ S ) 

+ (&£ 2 - 20£log(l + t/s) - ^tt 2 ^ Li 2 (-t/s) 

+ (8 1 + 20 log(l + t/s)) Li 3 {-t/s) 

+ 20 S 2 , 2 (~t/s) -201 Si, 2 (~t/s) - 28 Li 4 (-t/s) 

+ (28^-201og(l + t/ S ))C(3)-^}, (13) 
with t = \og(t/s). 

3. Differential equations for double-box 
master integrals 

The integrals I\ and Ii rr can now be chosen 
as a new basis of master integrals of the planar 
double-box topology, solving the problem of the 
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0(1/ e) coefficients of I\ and I2 in matrix element 
calculations [ |j . 

It is possible to construct a system of differen- 
tial equations satisfied by the double box master 
integrals [ ||, O. In terms of the new basis, these 
differential equations are [ || : 
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Expanding eqs. ( [14|Jl5| ) in e, and inserting the ex- 
pansion of Ii from ref. [0], of the pinched di- 
agrams (||) (see, eg. [ [l2]), and the result of 
cq. (p^), we find that they are indeed satisfied. 



Inspecting the right hand sides of the differen- 
tial equations, one notices that in eq. (|l4|), the 
coefficient of I] rr , and in eq. (|l^) , those of I\ and 
Iirr, are all proportional to d — 4. This means 
that, if h is known to O(e ), the O(e ) part of 
Ii rr is, a priori, only determined by the system of 
equations up to a ^-independent constant. 

One also observes that the system of differen- 
tial equations has a singular point at s + t = 0. 
This corresponds to the special kinematical con- 
figuration where p± + P3 — 0. At this point, the 
numerator of Ii rr becomes reducible: 

= (k +pi +p 2 +P3) 2 

= (k+P2? =Pi-P2+P 3 -s, (16) 



Pt 



so Iirr collapses to a linear combination of I\ and 
pinched diagrams. 

For a further discussion of this singular point, 
we refer to the contribution of Gehrmann and 
Remiddi [ who recently investigated the 

(equivalent) system of differential equations for 
I\ and 12- By excluding singular solutions of the 
differential equations as possible candidates, they 
were able to extract the order e piece of I\ — tl% 
using only the value of I\ at s + 1 = to O(e ) as 
input. The information contained in this quantity 
is equivalent to our eq. (13). 



4. 6 dimensions 

The master integrals I\ and I2 are both finite 
in d — 6 dimensions [ |lj| . This can be deduced 
from power counting considerations in momen- 
tum space; it is also easy to see by examining 
the arguments of the gamma functions in the 
Mellin-Barnes representation (pd|). Following the 
method of [ || , we relate these master integrals in 
d = 6 — 2e dimensions to master integrals of the 
new basis in d — 4 — 2e dimensions. Substituting 
e-expansions for the latter, we find that all pole 
terms indeed cancel, and the finite parts are 

rd=6 6 f fl i b 

1 i 



(s + t) t 
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ax = *^-6C(3), 
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(19) 
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a 2 = — + 6Li 4 (-t/s) + (I 2 + 7T 2 ) Li 2 (-t/s) 

+ 4£(C(3)-Li 3 (-t/ S )) + ^ 2 , (20) 

6 



b = f2C(3)-2Li3(-t/a)-y^)log(l + t/») 
+ l(^ +7r 2 )log2(1 + t/s) 



+ ^2£log(l + t/s)- — \ Li a (-t/«) 
+ 2ISi, 2 (-t/s)-2S 2 ,2(-V s )- 



(21) 



We checked these expressions by numerical in- 
tegration of Feynman parametric representations 
for the integrals in d = 6 in the region of negative 
s and t. This gives us an independent check on 
the finite parts of our result for I/ rr , eq. (^), and 
of Ji inref. [§. 

5. Conclusion 

We have resolved a technical issue related to 
the reduction of on-shell massless double box di- 
agrams to master integrals. Looking back, one 
can say that the basis of master integrals {Ii, I2] 
used in ref. [ || was an unfortunate choice for 
applications in d = 4 — 2e dimensions, because it 
leads to factors of 1/e in the reduction coefficients. 
The troublesome terms do not cancel in physical 
amplitudes [ ||. If one uses, instead, the basis 
{ii, Iirr}-, the factors of 1/e disappear. Our main 
result is the expansion of the integral Ij rr in terms 
of polylogarithms to O(e ), eq. (|l3|), which we 
derived from scratch using a Mellin-Barnes rep- 
resentation. (An equivalent result was obtained 
independently by Gehrmann and Remiddi [ |l3) , 
by solving a system of differential equations and 
using the first master integral from ref. [ |3) as in- 
put). Wc verified that our result, together with 
the first master integral, satisfies the differential 
equations ( fj"4|]T5| ) . As a further check on the con- 
stant which is not determined by the differential 
equations alone, we used our result to derive ana- 
lytical expressions for two finite double box in- 
tegrals in d = 6 dimensions, eqs. ( |l7| , |l8| ), and 
checked them by numerical integration. 
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